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Recent improvements in the experimental determination of properties of the Isovector Giant 
Quadrupole Resonance (IVGQR), as demonstrated in the A = 208 mass region, may be instrumental 
for characterizing the isovector channel of the effective nuclear interaction. We analyze properties 
of the IVGQR in 208 Pb, using both macroscopic and microscopic approaches. The microscopic 
method is based on families of non-relativistic and covariant Energy Density Functionals (EDF), 
characterized by a systematic variation of isoscalar and isovector properties of the corresponding 
nuclear matter equations of state. The macroscopic approach yields an explicit dependence of the 
nuclear symmetry energy at some subsaturation density, for instance S(p — 0.1 fm~ 3 ), or the neutron 
skin thickness Ar np of a heavy nucleus, on the excitation energies of isoscalar and isovector GQRs. 
Using available data it is found that S(p = 0.1 fm~ 3 ) = 23.3 ± 0.6 MeV. Results obtained with 
the microscopic framework confirm the correlation of the Ar np to the isoscalar and isovector GQR 
energies, as predicted by the macroscopic model. By exploiting this correlation together with the 
experimental values for the isoscalar and isovector GQR energies, we estimate Ar np = 0.14 ± 0.03 
fm for 208 Pb, and the slope parameter of the symmetry energy: L = 37 ± 18 MeV. 

PACS numbers: 21.60.Jz, 21.65.Ef, 24.30.Cz 



I. INTRODUCTION 

The Isoscalar Giant Quadrupole Resonance (ISGQR) 
was discovered in the 70s in inelastic electron and proton 
scattering experiments [IHH . For an experimental review 
we refer to [J|. Whereas the features of the low- lying 
quadrupole excitations depend on the number of parti- 
cles outside closed shells [fj - similarly to what occurs for 
the low-energy peak appearing in the isoscalar dipole re- 
sponse of neutron-rich nuclei [fj - the high-energy modes 
are expected to vary smoothly with the mass number A. 
In the case of the ISGQR, the excitation energy E^ can 
be estimated - considering the nucleus a quantal har- 
monic oscillator (QHO) - to be proportional to the shell 
energy-gap Hloq and, if the nuclear effective interaction 
is also velocity-dependent, to the nucleon effective mass, 
namely \Jmjm* (cf . Ref . Q ) • Because of this propor- 
tionality, the comparison of microscopic self-consistent 
calculations with experiments on the ISGQR has pro- 
vided valuable information on the value of m* Q , one of 
the most important quantities that characterize nucleons 
embedded in the nuclear medium 

At variance with the ISGQR, its isovector counterpart 
has remained elusive for quite a long time because of lack 
of selective experimental probes that can excite this res- 
onance. The accuracy in the experimental determination 
of the Isovector Giant Quadrupole Resonance (IV GQR) 
has been considerably improved only recently [10j . This 
important achievement will enable future measurements 
in different mass regions. The excitation energy of the 



IVGQR, E 1 /, is expected to vary smoothly with A. Op- 
posite to the ISGQR case, in the IVGQR neutrons and 
protons oscillate out of phase. Within the QHO assump- 
tion, the excitation energy of the high energy isovector 
mode should be correlated both with the shell gap (Hujo) 
and with the symmetry energy, as discussed below. Even 
though the symmetry energy S(p) is a basic component 
of the nuclear matter equation of state, it is still sig- 
nificantly undetermined [Ill4l4| |. At saturation density 
the symmetry energy is usually expressed in terms of its 
value, J = S(poo), and density slope, L = Sp 00 d p S(p)\ Poo . 
Also in the IVGQR case, for velocity dependent poten- 
tials parameterized in terms of an e ffective mass, the shell 
gap is modified as follows Huiq — > yJmJmFtiio^. 

In Section [TT] the theoretical basis of the non- 
relativistic Skyrmc and covariant Energy Density func- 
tionals (EDFs) is briefly presented. The formalism used 
in the present calculations is also outlined: mainly the 
Random Phase Approximation (RPA), and to some ex- 
tent the features of the Particle Vibration Coupling 
(PVC) approach. SectionQll]is divided into two parts. In 
Subsection IIII Al we analyze the strength functions and 
transition densities of the ISGQR and IVGQR in 208 Pb. 
In addition, the width of the IVGQR is evaluated using 
the PVC method. In Subsection IIII Bl we derive a macro- 
scopic model for the dynamics of the IVGQR. A detailed 
analysis of excitation energies of the ISGQR and IVGQR 
is performed employing two families of EDFs. Section ITVl 
summarizes the results and conclusions. 
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II. FORMALISM 



Mean-Field 



Self-consistent mean-field (SCMF) approaches to nu- 
clear structure have become increasingly complex and 
accurate. They represent an approximate realization of 
density functional theory (DFT) for atomic nuclei. This 
theory has been extensively applied to electronic sys- 
tems, based on the self-consistent Kohn-Sham scheme 
[15l - ll7l |. In nuclear physics different kinds of function- 
als are used, either local or non-local, based on a non- 
rclativistic (Skyrme, Gogny) or covariant representation. 
A review of modern SCMF models can be found in 
Ref . [ill . A common feature of these methods is their rel- 
atively simple structure. They are usually parameterized 
with about 10 constants adjusted to reproduce a selected 
set of ground state data. SCMF models yield accurate 
results for basic nuclear properties such as masses, radii 
and deformations, extending over the the entire chart of 
nuclides. 

Among the non-relativistic functionals, we employ here 
a set based on the Skyrme interaction [T|| [2(| ■ This in- 
teraction is of zero-range and density-dependent. One of 
the advantages of Skyrme functionals lies in the fact that 
the exchange terms (Fock terms) are simply proportional 
to the direct terms. 

Covariant SCMF, or rclativistic mean- field (RMF) 
models, have become another useful tool for the study 
of nuclear matter and finite nuclei [2ll - [2~4| ] . Nucleons are 
considered as Dirac particles coupled to effective mesons. 
The theory is Lorentz invariant and therefore preserves 
causality and provides a self-consistent description of the 
spin-orbit term of the nuclear effective force. Three ef- 
fective mesons comprise the minimal set necessary for a 
quantitative description of nuclear properties: a, oj, and 
p - in some cases the 5 meson has also been included (e.g. 
see Ref. [25|). Different types of effective Lagrangians 
have been considered. Well known examples are the 
Walecka-type models with linear and non-linear cr-meson 
self- interactions [2l[, such as the NL3 model [26| and, 
more recently, models based on density dependent finite- 
range meson- nucleon vertices, such as DD-ME2 [27|, or 
zero-range (point-coupling) interactions [Ifl [29| . 



B. Random Phase Approximations 

The introduction of a dynamical content into DFT- 
based models, leading to a time-dependent theory, is for- 
mally straightforward. In the realm of electronic density 
functionals, this scheme is called time-dependent density- 
functional theory (TDDFT) [30J ■ Nuclear physics imple- 
mentations exist, such as the time-dependent Hartree- 
Fock (TDHF) or time-dependent RMF. The lineariza- 
tion of the corresponding equations leads to the random 
phase approximation (RPA), in which collective nuclear 
excitations correspond to coherent superpositions of one 



particle-one hole (lp-lh) configurations. In particular, 
RPA is one of the most successful methods for the de- 
scription of nuclear excitations in the energy region of 
giant resonances (GRs). 

We briefly outline the basics of the discrete RPA for- 
malism [3l|,[32[. The RPA ground state is denoted |0), 
and \v) stands for a generic RPA excited state. For a 
given multipolc operator Fjm, the reduced transition 
probability is defined as 



B(EJ : 6-> i>) = \ (v\\Fj\\Q) 



P h 



. (1) 



where (p||i<j||/i) is the reduced matrix element of the op- 
erator Fjm, and X£ h and Y£ h are the RPA amplitudes. 
The strength function is defined by the relation 



S(E) = J2\{v\\Fj\\0) 



S(E - E v 



(2) 



where E v is the eigenenergy associated to the RPA eigen- 
state \v). The fc— moment of the strength function can 
be evaluated as follows: 



m k = J dEE k S(E) = J2\W\ P J 



16} 



(3) 



A useful quantity that provides information on the spa- 
tial features of the excited state is the transition density. 
Its integral with the radial part of a multipole operator 
yields the corresponding reduced transition amplitude for 
the given operator. For an RPA state \u) the radial part 
of the transition density, defined by 8p u (r) = {u\p(r)\0} = 
Spv(r)Yj M (f), is calculated using the expression: 

8p v (r) = ^=±=J2 
P h 



V2J+1 



i X ph 



Y ph. 



<p\\Yj\\h)^, (4) 



where u a (r) is the HF reduced radial wavefunction for 
the single-particle state a. The summation can run over 
proton and neutron states separately, thus defining the 
isoscalar (IS) and isovector (IV) transition densities: 

5pf(r) = 5pZ(r) + 5pP(r) 

8(F(r) = 5p2(r)-6f%(r). 

More details on our implementation of the non- 
relativistic and relativistic RPA models can be found in 

Refs. mm. 

The isoscalar and isovector quadrupole operators are 
defined by the following relations: 

A 



F2L = > 'riY 2 M(fi) 



E 



A 



= 1 



(G) 



(7) 
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C. Particle-vibration coupling 



III. RESULTS 



The SCMF approach to nuclear structure presents 
well-known limitations. For instance, it tends to un- 
derestimate the density of states around the Fermi en- 
ergy. Moreover, SCMF models cannot account for spec- 
troscopic factors of single-particle states, GR widths and 
decay properties. The nuclear field theory (NFT) (35l . [36j . 
based on the particle- vibration coupling (PVC), and in- 
troduced already in Ref. Q , provides a consistent frame- 
work for the treatment of beyond mean-field correlations. 
This framework allows for correlations between the static 
(single-particles) and the dynamic (phonons) parts of the 
mean-field. 

Recently a completely self-consistent approach to the 
PVC has been developed within the Skyrme framework 
[37j . The explicit expression for the basic vertex can be 
calculated using the RPA representation of the phonon 
state \n) with multipolarity L (equal to J as we deal only 
with natural parity phonons), thus obtaining: 



(i\\V\\j,nJ) = V2J+TJ2 X ph Vj(ihjp) 

ph (8) 
+ (-) jh - j " +J Y^ J Vj(ipjh), 

where Vj is the particle-hole J-couplcd matrix element 
Vj(ihjp) = Y.^-y^ mj+jh ~ mh ^ iTni ^ ~ mj\JM) 

{m} 



x (j P m p jh - m h \JM) 

x (jimij h m h \V\jjmjj p mp) . 



(9) 



For a detailed derivation of the reduced matrix element 
Eq. ©, we refer to the Appendix of Ref. j37j . 

In this work, in particular, we analyze the strength 
functions of the IVGQR. The coupling to low-lying vi- 
brations is the principal source of the GR width (the 
so-called spreading width) [38l - l40| . More information on 
the formalism that is used for the calculation of strength 
functions can be found in Refs. jU [42[. Here we just 
note the two main contributions to the spreading width: 
the self-energy of the particle (hole) that forms the reso- 
nance, i.e. the process in which a particle (hole) excites 
and reabsorbs a vibration, and the vertex correction that 
results from the exchange of a phonon between a particle 
and a hole (see below). 

Finally, it is important to note that the effective mass 
in the SCMF approach does not depend on the energy, 
and represents an average value for the whole nucleus as a 
function of the baryon density. It has been shown Q that 
taking into account the dynamical aspects of the nuclear 
mean-field is crucial for explaining the enhancement of 
the effective mass near the Fermi energy. 



GQRs: strength functions and transition 
densities 



In this subsection we analyze the main features of 
the strength functions associated with the isoscalar and 
isovector quadrupole response. Three Skyrme-type func- 
tionals, namely SAMi KDE [3, SkI3 and two 
relativistic functional, NL3 [26[ and DD-ME2 [27j] are 
considered. The Skyrme interactions have different ef- 
fective masses m*/m (0.68, 0.76, 0.58, respectively) and 
yield different values for the neutron skin thickness Ar np 



in 208 Pb (0.147 fm, 0.155 fm, 0.227 fm, respectively). The 
two covariant functionals are based on (i) finite-range 
meson exchange with non-linear self-interaction terms 
(NL3), and (ii) density-dependent meson- nucleon ver- 
tex functions (DD-ME2). Relativistic mean- field mod- 
els are known to yield rather low values for the non- 
rclativistic equivalent, or Schrodingcr effective mass, typ- 
ically around 0.6m at saturation density [HI, [46|. The 
NL3 functional predicts values of the neutron skin that 
are considerably larger compared to non-relativistic func- 
tionals, e.g., 0.279 fm for 208 Pb. The DD-ME2 functional 
yields the neutron skin thickness of 208 Pb: 0.193 fm. 

Figure Q] displays the isoscalar and isovector 
quadrupole transition strength functions, obtained by 
convoluting the RPA results with Lorentzian functions. 
The widths are taken in such a way that the total ex- 
perimental ISGQR and IVGQR widths are reproduced 
in the corresponding medium and high energy regions, 
respectively (see Table [TJ) . 

We start by analyzing the isoscalar quadrupole chan- 
nel in Fig. QJa). All models considered in the present 
study yield the ISGQR peak at excitation energies that 
are higher than the experimental value indicated by the 
arrow. It is well known Q that the energy of the isoscalar 
quadrupole resonance is closely related to the effective 
nucleon mass m* /m. Empirical ISGQR energies in heavy 
nuclei favor an effective mass close to 1. The effect of be- 
yond mean-field correlations on the isoscalar quadrupole 
strength functions were recently investigated in Ref. [42[ 
where, by using the SLy5 interaction |47j , a spreading 
width of the order of 2 MeV and a centroid energy of 
10.9 MeV were found, in very good agreement with data. 

The isovector spectrum shown in Fig. [TJb) consists of 
three distinct structures. The first one is the well known 
low-energy 2 + state at about 5 MeV, that we do not 
analyze in the present study. The second is the ISGQR 
that appears in the energy range between 10 and 15 MeV 
and, finally, the IVGQR located in the region above 20 
MeV. The two lower structures arise because of isospin 
mixing in the RPA states and, therefore, these could be 
excited both by isoscalar and isovector probes. In the 
high-energy region all interactions predict the existence 
of a collective IVGQR peak. Our results are in good 
agreement with experimental findings, both for the exci- 
tation energy of the IVGQR and the percentage of the 
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E (MeV) 

FIG. 1. (Color online) Isoscalar (a) and isovector (b) 
quadrupole strength functions. The strengths calculated with 
SAMi, KDE, SkI3, NL3 and DD-ME2 are shown. The exper- 
imental energies for the ISGQR (10.9 ± 0.1 MeV), and the 
IVGQR (22.7 ± 0.2 MeV) (weighted averages) listed in Table 
Uare indicated by arrows. 



energy-weighted sum rule (EWSR). The measured frac- 
tion for the latter is 56% [lOj . whereas theoretical pre- 
dictions range from 50% to 65%. Note that the EWSR 
fraction reported in Ref. [To[ refers to the classical ver- 
sion of the sum rule, that is, without the multiplicative 
factor (1 + kq), where kq is the isovector quadrupole 
enhancement factor [33| . 

More details about the structure of the IVGQR are 
provided by transition densities associated with the main 
peak of the isovector response. Fig.[5Ja) displays the neu- 
tron and proton transition densities, and in Fig. [2Jb) we 
plot the corresponding isoscalar and isovector transition 
densities calculated with the functionals SAMi and DD- 
ME2. The other functionals considered in this work yield 
similar transition densities and we do not show their re- 
sults. The positions of the proton (r p ) and neutron (r n ) 
root mean square (rms) radii correspond to the edges of 
the shaded region that, in this way, denotes the neutron 
skin thickness calculated with a given functional. For all 
functionals and, in particular for those used in Fig.[2j one 
notices that protons and neutrons yield similar contribu- 
tions but with opposite signs to the transition densities 
in the surface region. This shows that the excitation is 



predominantly isovector. In the bulk of the nucleus one 
finds a non negligible isoscalar component, even when the 
state is mainly isovector. This is, in particular, the case 
for the relativistic DD-ME2 functional. 

Although the aim of this work is the study of global 
properties, and in particular the energy of the IVGQR 
using relativistic and non-rclativistic EDFs, we com- 
plete the theoretical analysis by calculating the width 
of this important resonance. The model is described in 
Section III Gl and takes into account beyond mean-field 
correlations. In Fig. [3{a), we display in a diagram- 
matic way the processes in which a particle (hole) excites 
and reabsorbs a vibration (left), and the vertex correc- 
tion in which a particle and a hole exchange a phonon 
(right). In Fig.^b) the probability of finding the isovec- 
tor quadrupole resonance state per unit energy is plotted, 
calculated with the SAMi functional. Phonons with mul- 
tipolarity L = 0,1,2,3 and natural parity are included 
in the model space. The RPA model space in this cal- 
culations is taken as in Ref. [42j |. With this choice, the 
EWSRs are satisfied up to about the 99%. As in Ref. [Hj], 
we impose a lower cutoff on the collectivity of the inter- 
mediate RPA states for two reasons: firstly the RPA docs 
not provide a good description of non-collective states 
and, secondly, phonons with low collectivity would ne- 
cessitate taking into account corrections that arise from 
the Pauli principle. For these reasons only phonons with 
energy lower than 30 MeV and a fraction of EWSR larger 
than 5% are included as intermediate states. The sin- 
gle RPA state splits into two components: the peak at 
higher energy is barely affected by increasing the num- 
ber of intermediate phonons, whereas the one at lower 
energy broadens and is shifted downwards as the number 
of phonons increases. The energy centroid is shifted by 
1 MeV, from 22 MeV for the RPA to 21 MeV for the 
PVC model. The spreading width is 3.8 MeV. Both re- 
sults are in reasonable agreement with the experimental 
results listed in Table Q] 

B. Sensitivity of the GQRs to the symmetry 
energy and the nucleon effective mass 

Based on the QHO approach [tJ , in the Introduction we 
have discussed the main characteristics of excitation ener- 
gies of isoscalar and isovector GQRs. In the ISGQR case 
the excitation energy can be estimated from the QHO 
formula, 

Ef = y/2m/m*hu; , (10) 

where the last factor is the shell gap ~ 41A -1 / 3 . Such 
a semi-empirical approach reproduces the average trend 
and predicts the correlation between the ISGQR ex- 
citation energy and the nucleon effective mass that is 
also confirmed by microscopic calculations Q. For the 
IVGQR the QHO is expected to provide a direct relation 
between the excitation energy and the isovector proper- 
ties of a nuclear effective interaction. 
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FIG. 2. (Color online) Neutron and proton (a) and isoscalar and isovector (b) transition densities for the main peak of the 
isovector response, as a function of the radial distance. The predictions for the SAMi and DD-ME2 functionals are shown. The 
proton (r p ) and neutron (r„) rms radii are indicated by the edges of the shaded region. 



To show this relation we derive a macroscopic for- 
mula for the excitation energy of the IVGQR that ex- 
plicitly exhibits the connection to both the effective mass 
and the symmetry energy at some subsaturation density 
S(p = 0.1 fm -3 ). In a first step we bypass the depen- 
dence on the effective mass and replace it with the IS- 
GQR energy. In this way the value for S(p = 0.1 fm~ 3 ) 
will be determined from experimental results only. Sub- 
sequently a quantitative investigation of this correlation 
is performed for the case of 208 Pb by employing families 
of EDFs. 

We use available data on ISGQR and IVGQR in the 
A = 208 mass region (Table HJ to estimate the value of 
the symmetry energy at 0.1 fm _1 . Although the IVGQR 
in Ref. [l(| was measured in 209 Bi, calculations are car- 
ried out for 208 Pb. The difference in energy of the nuclear 
response of 209 Bi and 208 Pb should scale with A~ 1/3 0, 
that is, it should be smaller than a few %o. Another im- 
portant reason for limiting the study to 208 Pb is that it 
is a spherical double magic nucleus, and thus the depen- 
dence on the effective mass or the symmetry energy will 
not be screened by deformation or paring effects. 

The EDFs we employ in this study are based on differ- 
ent theoretical frameworks. One is the non-relativistic 
Skyrme-Hartree-Fock approach (SAMi HU), and the 
other is the relativistic mean-field with density depen- 
dent meson- nucleon vertices (DD-ME [55j). We have con- 
sidered families of functionals with systematically varied 
properties in the isoscalar and isovector channels. For 
SAMi, using the fitting protocol described in the origi- 
nal reference j4[|, we have first fixed the values of the 
nuclear incompressibility (iVx> = 245 MeV) and the ef- 
fective mass (m*/m = 0.675), whereas the values of the 
symmetry energy at saturation (J) have been varied from 
27 MeV (SAMi-J27) to 31 MeV (SAM-J31) in steps of 
1 MeV. Then, by fixing the values of K x = 245 MeV, 
J = 28 MeV and L — 44 MeV, we have varied the effec- 



TABLE I. Data for the IVGQR and ISGQR in 208 Pb. (a) 
These experimental values are for 209 Bi, and the EWSR corre- 
sponds to the classical value (see text). Weighted averages 
obtained from the data. 
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IVGQR 


24.3±0.4 


4.5±0.5 


140 


[48] 




22.5 


9 


100 


[49] 




20.2±0.5 


5.5±0.5 


140±30 


[50] 




23.0±0.2 


3.9±0.9 


56± 6 
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Average 


22.7±0.2 


4.8±0.3 
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ISGQR 


10.60±0.25 


2.8±0.25 


100 


[51] 




11 ±0.2 


2.7±0.3 


105±25 


[52] 




10.9 ±0.3 


3.1±0.3 


120-170 
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11.0 ±0.3 


3.3±0.3 


100-150 


[53] 




10.9 ±0.3 


3.0±0.3 


100±13 


[5J 


Weighted 








Average 


10.9 ±0.1 


3.0±0.1 







tive mass from m*/m = 0.65 (SAMi-m65) to 0.85 (SAMi- 
m85) in steps of 0.05. In the case of the relativistic func- 
tionals DD-ME, we have adopted the set of interactions 
introduced in Ref. [55j], in which J was systematically 
varied from 30 MeV to 38 MeV in steps of 2 MeV (sets 
from DD-MEa to DD-MEc). This kind of analysis allows 
to identify possible correlations. 



1. Macroscopic model for the excitation energy of the 
Isovector Giant Quadrupole Resonance 

In the QHO model (see Eqs. (6-379) and (6-381) 
in Ref. 0) the excitation energy of the isovector giant 
quadrupole excitation mode can be written in the follow- 
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with the parameters characterizing the nuclear symmetry 
energy S(p) can be found, within the SCMF approach, 
by using the empirical law of Ref . [HJ , where it has been 
demonstrated that the symmetry energy of a finite nu- 
cleus af^(A) equals the symmetry energy S(p) of the 
infinite system at some sub-saturation density pa - ap- 
proximately 0.1 fm~ 3 for the case of heavy nuclei such as 
208 Pb. Hence, one can rewrite Eq. (fTTj) as 
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FIG. 3. (Color online) Diagrams contributing to the strength 
function of the GR [panel (a)]. Probability P to find the 
IVGQR state at an energy E [panel (b)]. Different curves 
are obtained when the phonons listed in the legend are used 
as intermediate states. The label RPA [black-dashed line in 
panel (b)] refers to the curve calculated in the RPA with a 
Lorentzian width of 1 MeV. 



ing form 
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IV 



5 H 2 V sym (r 2 ) 
4 2m (hajofir 4 )' 



(11) 



where V sym is the symmetry potential proportional to the 
liquid drop model (LDM) parameter b^°^: V sym — 4b p ° I ^ 
[see Eq. (2-28) in Ref. 0]. &P y * n can be written as bf°^ = 

frsym — &sym ~ ^sym — 2S' kln (p 00 ). In the non-relativistic 
approximation the kinetic contribution to the symmetry 
energy at nuclear saturation is S kln (p aD ) « ef^/3 [see 
Eq. (2-13) in Ref. 0], where e Foo = h 2 klj2m ~ 37 
MeV is the Fermi energy for symmetric nuclear matter 
at saturation density. The relation that connects 6 sym 
with the "standard" liquid drop parameter a^^ 1 reads 
b sym « 2a^VM ( see Eq (2-12) from Ref. 0). Since giant 
resonances in finite nuclei are not pure volume modes, 
it is important to take into account surface corrections 
and, therefore, one may identify a]^ 1 with the Droplet 
Model (DM) parameter that contains surface corrections 
(A) [56|, |57| ■ The connection of the latter quantity 



,DM 
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IV 



|(7ttJ ) 



2 , r £ F 

6 



6 



A 2/3 

A 2 / 3 L 



[S(pa) 
S(pa) 



gkin 



(Poo)]} 



1/2 



1/2 



(12) 



where we have approximated the factor ^ {^) 2 ^ = 
2.0113 by 2 on the r.h.s., and considered (r n ) = 
3r%A n / 3 /{n + 3) where r = [3/(4tt / 9 00 )] 1 / 3 . As in the 
case of the ISGQR, for velocity dependent potentials pa- 
rameterized in terms of an effective mass, the shell gap 
is modified as follows 



E 



IV 



m \ 2 

TO* 



J A 2 ' 3 \ £ Foo 
or, equivalently, by using Eq. (|10[) 



1/2 



(13) 



E 



IV 



A 2 / 3 \ e Fo 



- 1 



~\ 1/2 



(14) 



From Eq. (|13|) and approximating fiajQ 
note some interesting features: 



41A" 1 /3 i we 



sym 



• E~y depends, as its isoscalar counterpart, on the 
effective mass at saturation and, in addition, on 
the symmetry energy at some sub-saturation den- 
sity pa (the Fermi energy at saturation can be con- 
sidered as constant compared with the variation of 
other quantities) . E 1 ^ increases for decreasing val- 
ues of m*, and increasing values of S(pa)- 

• The larger the neutron skin thickness in a heavy 
nucleus such as 208 Pb, the lower the excitation en- 
ergy of the IVGQR. This characteristic can be un- 
derstood as follows. If one expands S(p) around the 
nuclear saturation density as S(p) w J — Le, where 
e = (poo — p)/p, it can explicitly be shown that at 
the sub-saturation density pa, fixing £7* to the ex- 
perimental value and for small variations of J, E^ 
decreases for increasing values of L. The latter is 
linearly correlated with the neutron skin thickness 
(llTtl4| . which increases with L (see below). 

One of the most important consequences of our ap- 
proach is that from Eq. (TT4")) one can find a relation that 
expresses S(p A ) in terms of E™ , E^, and the Fermi en- 



7 



ergy at nuclear saturation ep 

A 2/3 



that is, 



S( P/ 



24 4 



£ F 



(^) -2(^) J+^oc) 
A 2/3 



2(El s Y +1 .(15) 



By inserting the weighted averages of the experimental 
values for = 22.7 ± 0.2 MeV and £^ s = 10.9 ± 0.1 
MeV (see Table [J), and by using pa=208 = 0.1 fm~ 3 , we 
find 5(0.1 fm. -3 ) = 23.3 ± 0.6 MeV, in very good agree- 
ment with the estimate reported in Ref. J58|: 23.3 MeV 
< 5(0.1 fm~ 3 ) < 24.9 MeV. Note that the quoted error 
does not include an estimate of the theoretical uncer- 
tainty. 

Since we are also interested in determining correlations 
between observable quantities, we elaborate on Eq. ([15)) 
in order to explicitly relate the excitation energies of the 
isoscalar and isovector GQRs with the neutron skin thick- 
ness of a heavy nucleus. For that, we use the DM expres- 
sion for the neutron skin thickness that can be written 
as follows UM 



A, 



np 



A surf 
L - x ' np 



( r 2)l/2 



1 



J 



(16) 



where I = (N — Z)/A is the relative neutron excess, 
I c = e 2 Z/(20JR) and Ar™ rf is the surface contribution 
to the neutron skin thicknes^]. The latter, for the case 
of 208 Pb, has a value of 0.09 fm [Hi, if calculated with 
a large set of EDFs. Combining Eqs. p^|) and (|16l) one 
finds, 



Ar: 



surf 



(r 



2\l/2 



3 Ic 



A 2 / 3 



71 

2 



J 



.(17) 



This expression explicitly relates the neutron skin thick- 
ness of a heavy nucleus with the corresponding GQRs 
energies, and these can directly be determined in exper- 
iment. Within our approach, only the parameter J and 
Ar^ rf contain a non-negligible uncertainty. The appro- 
priate value of J to be used in the expression above can 
be deduced from the systematic analysis carried out in 
Refs. [3, |63|. The weighted average of the constraints 
considered in [3 yields J = 32.4 ± 0.4 MeV, but in the 
following we adopt a somewhat larger uncertainty J = 
32 ± 1 MeV. For the case of 208 Pb, Ar™ rf = 0.09 ± 0.01 
fm is consistent with the microscopic calculations of 
Ref. Using Eq. dTTJ) and the data for the GQRs 

energies, we find Ar„ p = 0.22 ± 0.02 fm. This value is 
close to the upper limit derived from available estimates 
Ar np = 0.18 ±0.03 fm 
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FIG. 4. (Color online) Excitation energy of the ISGQR in 
208 Pb as a function of i/m/m*, calculated with the SAMi-m 
and SAMi-J family of functionals. On the horizontal upper 
axis we display the corresponding values for m* /m. The data 
from Table U are also included (square and shaded band). 



2. The SAMi and DD-ME families of functionals 

In Fig.@]we display the exc itation energy of the ISGQR 
in 208 Pb as a function of y^m/m*, calculated with the 
SAMi-m and SAMi-J families of functionals. The plot 
nicely illustrates the well known correlation between 
and yjm/m* . It also shows that the variation of E 1 ^ for 
the SAMi-J family - for which J ranges from 27 to 31 
MeV, with a fixed value of m* , is smalo 

Fig. [5] displays the predictions of the SAMi-m and 
SAMi-J families of functionals for the excitation energy of 
the IVGQR in 208 Pb as a function of the effective mass 
(Fig. (5][a)) and neutron skin thickness (Fig. [5jb)), and 
those of the DD-ME family for the excitation energy of 
the IVGQR as a function of the neutron skin thickness 
(Fig. EJc)). Based on the discussion above, we expect a 
linear correlation between (S^X) and cither m/m* or 
Ar„ p . These two correlations are clearly displayed by 
the linear fits in Figs. [SJa) and EJb) for the SAMi func- 
tionals. In Fig. [5jc), consistent with the macroscopic 
model of Subsection IIII B II and with the non-relativistic 
results, a linear correlation is found between the square 
of the excitation energy of the IVGQR in 208 Pb and the 
neutron skin thickness predicted by the DD-ME family 
of functionals. 

These results, as well as the macroscopic model of Sub- 
section IIII B 11 show that a measurement of the excita- 
tion energy of the IVGQR in 208 Pb determines only a 
combination of the excitation energy of the ISGQR, [see 
Eqs. (TT3"]) and (TT4")) ] and Ar np (or the value of the slope 



1 Ar^ rf = vW5 [5 ( b l - ftp) / (2K)] , where b n and b p are the 
surface widths of the neutron and proton density profiles, re- 2 The neutron radius increases with J |61| . and a larger size of 
spcctively. the nucleus implies a lower ISGQR excitation energy. 
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FIG. 5. (Color online) Square of the excitation energy of the 
IVGQR in 208 Pb as a function of the effective mass (a), and 
the neutron skin thickness (b,c), predicted by the SAMi-m 
and SAMi-J (panels (a) and (b)), and DD-ME (c) families of 
energy density functionals. 



of the symmetry energy at saturation L, see Sec. IHIB 11 
and cf. Refs. |llT - fl4| ). but not their individual values. 

In Table Hi] we also compare the values for the sym- 
metry energy at 0.1 fm -3 S(pa) obtained in asymmetric 
nuclear matter with the SAMi-J and DD-ME families of 
EDFs, to the corresponding values S(pa) calculated us- 
ing Eq. (fl~5j) . As explained in subsection MI B 1[ the ex- 
pression for S(pa) is based on a (non-relativistic) QHO 
approximation and depends on the excitation energies of 
the GQRs and the kinetic contribution to the nuclear 
symmetry energy at saturation. We find that the differ- 
ence - calculated as the rms deviation - between S(pa) 
and S(pa) is less than 1 McV for the SAMi-J family, 
whereas it is around 3 MeV for the DD-ME functionals. 

Eq. (HH) shows that [(E™) 2 - 2{Eff]/J is linearly 
correlated with Ar np . This correlation is illustrated in 
Fig. [S] where we plot the Ar np calculated with the SAMi- 
J and DD-ME functionals, as functions of {(E™) 2 — 
2{E l x) 2 ]/ J. Both families show a high linear correla- 
tion (r = 0.98) between these two quantities, but predict 
different slopes. The slope obtained in the macroscopic 
model is: (r 2 )^ 2 (I - 7 c )A 2 / 3 /(36£ F J [cf. Eq. fTQ], 
independent of 5 kin ( / 9 oc ). For 208 Pb this yields 0.025 
MeV -1 fm, in very good agreement with the value 0.027 
McV -1 fm found for the SAMi family. The macroscopic 
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FIG. 6. (Color online) Values of Ar np in 208 Pb as functions 
of {(E™) 2 - 2{Eff]/J, calculated with the SAMi-J and DD- 
ME functionals. The dashed line and shaded band indicate 
the experimental value and corresponding uncertainty (see 
text). 



formula obviously does not apply to the relativistic case 
since the slope for the DD-ME family of functionals is 
0.057 MeV- 1 fm. 

Using the linear correlations shown in Fig. |6l the ex- 
perimental values for E 1 / = 22.7 ± 0.2 MeV and = 
10.9 ±0.1 MeV from Table H and the value J = 32 ± 1 
MeV that yields {(E™) 2 - 2{Ef ) 2 ]/ J = 8.7 ± 0.4 MeV, 
one finds Ar np = 0.14 ± 0.03 fm for the DDME fam- 
ily of functionals, and from the analysis of the SAMi- 
J functionals Ar np = 0.14 ± 0.02 fm. The total range 
of allowed values 0.11 fm < Ar np < 0.17 fm is rather 
broad but in reasonable agreement with previous stud- 



ies: Ar np = 0.18±0.03fm [Ij, and Ar np = 0.188±0.014 
fm [62j . Finally, this result for the neutron skin thickness 
of 208 Pb allows us to estimate the value of the slope of 
the symmetry energy at saturation for the DDME and 
SAMi-J families. Fig. [7] shows that this value is in the 
interval: 19 McV < L < 55 MeV. Wc note that the cor- 
relation coefficient is rather hig h and in agreement with 
those obtained in Refs. [IIIptI ]. Our constraint on L is 
also in a gree ment with previous estimates (see Fig.l in 

Refs. [T3J13, H3l)- 



IV. SUMMARY AND CONCLUSIONS 

Motivated by recent experimental developments fl"ot | . 
in the first part of this work we have investigated the 
excitation energy of the IVGQR in 208 Pb using a set of 
non-relativistic and covariant EDFs. The theoretical re- 
sults are in good agreement with the experimental find- 
ings for the excitation energy and the EWSR. We have 
also analyzed the transition densities associated to the 
principal RPA states that correspond to this resonance. 
It has been found that the isovector character is dominant 
in the surface region of the nucleus, whereas the interior 
part displays a non-negligible isoscalar component. The 
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TABLE II. Theoretical values for the neutron skin thickness Ar np of Pb, the symmetry energy S(p) at p = 0.1 fm and at 
poo (J), the kinetic contribution to the symmetry energy, S kln (poo J3, and the Isoscalar and Isovector GQR excitation energies 
E^f and , respectively. We also display 5(0.1 fm~ 
of fm and all other quantities are in units of MeV. 



as predicted by Eq. (|15[1 [5*(0.1 fm ' )]. The values of Ar np are in units 
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FIG. 7. (Color online) Neutron skin thickness Ar np of 208 Pb 
as a function of the slope parameter of the symmetry energy at 
saturation density, for the two families of functionals: S AMi- J 
and DDME. 



deduce, from the measured excitation energies of the IS- 
GQR and IVGQR, the symmetry energy at a subsatu- 
ration density 0.1 fm -3 (or the neutron skin thickness 
Ar np of the considered nucleus). The estimated value 
5(0.1 fm -3 ) = 23.3 ± 0.6 MeV, is in very good agree- 
ment with previous findings [HH . A strong correlation be- 
tween (E^y - 2 (^ s ) j J and Ar np has been found 

for the two families of EDFs considered in this work. This 
means that data on the excitation energy of the ISGQR 
and the IVGQR can be used to determine the neutron 
skin thickness of a heavy nucleus and the slope of the 
symmetry energy at saturation. With this approach we 
have obtained for the neutron skin thickness of 208 Pb: 
Ar np = 0.14 ± 0.03 fm, and the slope parameter of the 
symmetry energy L = 37 ± 18 MeV. These values are 
compatible with previous estimates. 



spreading width of the resonance has been calculated us- 
ing the Particle Vibration Coupling approach, and the 
resulting value is in good agreement with the data. 

In the second part we have focused on the relation be- 
tween the excitation energy of the GQRs and the isovec- 
tor properties of effective nuclear interactions. For this 
purpose a macroscopic formula has been derived, based 
on the quantal harmonic oscillator model Q and the ap- 
proach of Ref. fl2j |. Despite its simplicity, this formula 
provides a connection between the macroscopic picture 
of the IVGQR and microscopic calculations based on ac- 
curately calibrated families of EDFs. 

Using the analytic expression we have been able to 



ACKNOWLEDGMENTS 

The authors would like to thank M. Centelles for use- 
ful discussions. We acknowledge support of the Italian 
Research Project "Many-body theory of nuclear systems 
and implications on the physics of neutron stars" (PRIN 
2008); MZOS - project 1191005-1010 and the Croatian 
Science Foundation; Grants Nos 10875150 and 11175216 
of the National Natural Science Foundation of China and 
the Project of Knowledge Innovation Program (PKIP) of 
Chinese Academy of Sciences, Grant No. KJCX2-EW- 
N01. 



10 



R. 



Pitthan 



and 



T. 



|Phys. Lett. B 36, 563 (1971) 



Fukuda 



and 



|Phys. Rev. Lett. 29, 1109 (1972) 



M. Lewis and 

|Nucl. Phys. A 196, 337 (1972) 



Walcher, 



Torizuka, 



Bertrand, 



M. Harakeh and A. van der Woude, 
Giant Resonances. Fundamental High-Frequency Modes of 



(Clarendon Press, Oxford, 2011). 
N. Tsoneva and 



H. 



Lenske, 



Physics Letters B 695, 174 (2011) 



X. Roca- Maza, G. Pozzi, M. Br enna, K. Mizuyama, and 

G. Colo, |Phys. Rev. C 85, 1 (2012)] 

A. Bohr and B. R. Mottelson, Nuclear Stucture , Vol. I & 

II (W. A.Ben jamin Inc., 1975). 

J.-P. Blaizot, |Phys. Rep. 64, 171 (1980)] 

C. Mahaux, P. F. Bortignon, R. A. Broglia, and C. H. 

Dasso, Phys. Rep. 4, 1 (1985). 

S. S. Henshaw, M. W. Ahmed, G. Feld- 
man, A. M. Nathan, and H. R. Weller, 
|Phys. Rev. Lett 107, 222501 (2011) 



B. Alex Brown, |Phys. Rev. Lett. 85, 5296 (2000)| 
M. Centelles, X. Roca-Maza, X. Vinas, and M. Warda, 
|Phys. Rev. Lett. 102, 122502 (2009)] 
X. Roca-Maza, M. Centelles, X. Vinas, and M. Warda, 
|Phys. Rev. Lett. 106 252501 (2011)' 



M. B. Tsang, J. R. Stone, F. Camera, P. Danielewicz, 
S. Gandolfi, K. Hebeler, C. J. Horowitz, J. Lee, 
W. G. Lynch, Z. Kohley, R. Lemmon, P. Moller, 
T. Murakami, S. Riordan, X. Roca-Maza, F. Sammar- 
ruca, A. W. Steiner, I. Vidana, and S. J. Yennello, 



|Phys. Rev. C 86, 015803 (2012) | 

M. Brack, Helv. Phys. Acta 58, 715 (1985). 

R. N. Schmid, E. Engel, and R. M. Dreizler, 

|Phys. Rev. C 52, 164 (1995)] 

R. N. Schmid, E. Engel, and R. M. Dreizler, 



Phys. Rev. C 52, 2804 (1995) 



M. Bender, P.-H. Heenen, 
|Rev. Mod. Phys. 75, 121 (2003) | 
D. Vautherin and D 
|Phys. Rev. C 5, 626 (1972)] 
M. Beiner, H. Flocard, N. 



and P.-G. Reinhard, 



M. 



Brink, 



V. Giai, and P. Quentin, 



|Nucl. Phys. A 238, 29 (1975)| 
B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 
(1986). 

P.-G. Reinh ard, |Rep. Progr. Phys. 52, 4 39 (1989) 
B. D. Serot, [Rep" Progr. Phys. 55, 1855 (1992)| 
D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, 
P. Ring, |Phys. Rep. 409, 101 (2005)| 
X. Roca-Maza, X. Vinas, M. Centelles, P. Ring, 
P. Schuck, [Phys. Rev. C 84, 054309 (2011)] 
G. A. Lalazissis, J. Konig, and P 
|Phys. Rev. C 55, 540 (1997)] 
G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring 



and 



and 



Ring, 



[31] D. J. Rowe, Nuclear Collective Motion (Methuen and Co. 
Ltd., 1970). 

[32] P. Ring and P. Shuck, The Nuclear Many-Body Problem 

(Springer, 2004). 
[33] G. Colo, L. Cao, N. V. Giai, and L. Cape lli, 

[Computer Physics Communications 184, 142 (2013) | 
[34] N. Paar, D. Vretenar, E. Khan, and G. Colo, Rep. Prog. 
NuclefflpsLcIBafrifegLS (2007). 

D. R. Bes, G. G. Dussel, R. A. Broglia, R. Liotta, and 
B. R. Mottelson, |Physics Letters B 52, 253 (1974)| 



[35 
[36 
[37 
[38 
[39 

[40 

[41 
[42 
[43 
[44 
[45 
[46 

I 47 
[48 
[49 

[50 
[51 

[52 

[53 
[54 



|Phys. Rev. C 71, 024312 (2005) 




[55] 


W.-H. Long, N. V. Giai, and J. 


Meng, 




|Physics Letters B 640, 150 (2006) 




[56] 


T. Niksic, D. Vretenar, and P. 


Ring, 




|Phys. Rev. C 78, 034318 (2008) | 




[57] 


E. Runge and E. K. U. 


Gross, 


[58] 


|Phys. Rev. Lett. 52, 997 (1984) | 







P. F. Bortignon, R. A. Broglia, D. R. Bes, and R. Liotta, 
|Phys. Rep. 30, 305 (1977)] 

G. Colo, H. Sagawa, and P. F. Bortignon, 



Phys. Rev. C 82, 064307 (2010) 



R. 



A. 



Brc 



P. F. Bortignon and 
|Nucl. Phys. A 371, 405 (1981)| 

P. F. Bortignon^ ~K~. Bracco, and R. A. Broglia, 
Giant Resonances: Nuclear structure at finite temperature , 
Contemporary Concepts in Physics (Harwood Academic 
publisher, 1998). 

G. F. Bertsch, P. F. Bortignon, and R. A. Broglia, 
|Reviews of Modern Physics 55, 287 (1983)| 
G. Colo, N. Van Giai, and P. F. 
[Physical Review C 50, 1496 (1994) [ 
M. Brenna, G. Colo, and 
[Phys 



P. F. 



f s. Rev. C 85, 014305 (2012) 

^G^ 



X. Roca-Maza, 



Colo, 



and H 



Bortignon, 
Bortignon, 
Sagawa, 



|Phys. Rev. C 86, 031306(R) (2012) [ 
TJ! K. Agrawal, 5! Shlomo, 
|Phys. Rev. C 72, 014310 (2005) | 
P.-G. Reinhard and 

|Nucl. Phys. A 584, 467 (1995) 



and V. K. Au, 



H. 



M. Jaminon and 

|Phys. Rev. C 40, 354 (1989) | 
E. Chabanat, P. Bonche, P. Haensel, J. 



Flocard, 



Mahaux, 



Meyer, and 



R. Schaeffer, |Nucl. Phys. A 635, 231 (1998) [ 

R. Leicht, M. Hammen, K. Schelhaas, and B. Ziegler, 

|Nucl Phys. A 362, 111 (1981)[ 

K. Schelhaas, T. Henneberg, M. Sanzone-Arenhvel, 
N. Wieloch-Laufenberg, U. Zurmhl, B. Ziegler, M. Schu- 
macher, and F. Wolf, |Nucl. Phys. A 489, 189 (1988)] 
D. Dale, R. Laszewski, and R Alarcon, 



Phys. Rev. Lett. 68, 3507 (1 992) 



M. Buenerd, in Journal de Physique Collogue , Vol. 45 
(France, 1984) pp. 115-34, international Symposium on 
Highly Excited States and Nuclear Structure, 5-8 Sept. 
1983, Orsay, France. 

D. H. Youngblood, P. Bogucki, J. D. Bron- 
son, U. Garg, Y. W. Lui, and C. M. Rozsa, 



Phys. Rev. C 23, 1997 (1981)] 



S. Brandenburg, Ph.D. thesis, University of Groningcn 
(1985). 

D. H. Youngblood, Y.-W. Lui, H. L. Clark, B. John, 



Y. Tokimoto, and X. Chen, |Phys: Rev. C 69, 1 (2004) 
D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C 68, 
024310 (2003). 

W. D. Myers and W. Swiatecki, 

|Ann. Phys. 55, 395 (19"69)1 

W. Myers and W. Swiatecki, |Ann. Phys. 84, 186 (1974) 
L. Trippa, G. Colo, 
|Phys. Rev. C 77, 061304 (2008) | 



and 



Vigezzi, 



11 



[59] M. Centelles, X. Roca-Maza, X Vinas, and M. Warda, [62] B. K. Agrawal, J. N. De, and S. K. Samaddar, Phys. 
|Phys. Rev. C 82, 054314 (2010) | Rev. Lett, in press (2012). 

[60] J. M. Lattimer and Y. Lim, [63] X. Vinas, M. Centelles, X. Roca-Maza, and M. Warda, 

| "Constraining the symmetry parameters of the nuclear interactic|iA,rF| Conf. Proc. pp. -104; doi:http://dx.doi.org/ 1491, 101 (2 012) 



(2012), |arXiv7nucl-th/1203.4286 



[61] R. Furnstahl, |Nuclear Physics A 706, 85 (2 002) 



